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A discrete spacetime is a random partially ordered set. It is obtained by sam- 
pling n points from a compact domain in a spacetime manifold w.r.t. the volume 
form. The partial order is given by restricting the causal structure on the domain 
to the sampled points. These objects were introduced by physicists in an attempt 
to reconcile gravity with quantum mechanics, [401 1451 119] . They received very 
little direct attention in the mathematics literature. 

Recently, discrete spacetimes have appeared in very diverse applications which 
are very far removed from their original motivation. In this self-contained survey we 
will present some results and numerical calculations regarding the asymptotic be- 
havior of these objects. Some of the more basic results are stated as theorems since 
they are particularly useful in the applications. Numerical calculations are pre- 
sented with the hope that they will incite some theoretical work, aimed at proving 
some related assertions. We then present a few of the applications by shamelessly 
following the pattern of the recent survey paper [21j , which describes some related 
material. We will present a short list of seemingly unrelated problems and show 
how they fit into the discrete spacetime formalism. While this is a survey paper, it 
does contain some new material, namely, the statement of theorem ll.31 the (simple) 
solution to the dimension reconstruction problem, the numerical calculations of fig- 
ure [21 the connect-the-dots application and the relation between discrete spacetime 
and the literature on maximal layers. A more comprehensive and extensive treat- 
ment of the subject along with full proofs of the basic results will be presented in 
a forthcoming publication. 

Our own research relating to discrete spacetime owes much to the constant 
encouragement which was provided by Percy Deift and it is our sincere pleasure to 
dedicate this paper to him on the occasion of his 60th birthday. 



1. The asymptotic behavior of discrete spacetime 

A spacetime (Lorentzian manifold) Af is an to = c? -f 1 dimensional real, con- 
nected, C°° manifold with a globally defined C°° symmetric tensor field g, of type 
(0,2), which is non-degenerate and Lorentzian. By Lorentzian we mean that for 
any a ^ M there is a basis of Ta(M), the tangent space to M at a relative to 
which ga is given by the diagonal matrix with diagonal entries (—1,1,1,...,!). We 
consider a set of local coordinates xq, . . . ,Xd on the manifold. In terms of the local 
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coordinates it is convenient to consider the expression 
(1) ds^ — gi^j dxi dxj 

i j 

The metric defines a canonical volume form which, using the coordinates, may be 
written as 



(2) Vol, = V|det(5a)|d™x 

The metric ga also classifies the tangent vectors v G TaM into timelike vectors, 
if ga{v,v) > 0, null or lightlike, if ga{v,v) — and spacelike ii ga{v,v) < 0. The 
set of timelike vectors consists of a double cone. We say that a spacetime is time 
orientable if there is a non-vanishing vector field V consisting of timelike vectors. 
We will assume that all our spacetime manifolds are time orientable. An orientation 
provides a continuous choice of one out of the two cones (the one containing the 
vector Va) and we declare this cone to be the future-pointing cone. A parametrized 
curve in the manifold is said to be future-pointing timelike if all its tangent vectors 
belong to the future-pointing cone. 

The simplest spacetime manifold is flat Minkowski space, which is R™ equipped 
with the constant metric 

(3) ds^ = dxQ — dx\ — dx\ — ... — dx'\ 

Other natural examples to consider are the constant curvature hyperbolic (anti de- 
Sitter) and spherical (de-Sitter) spacetime manifolds. In two dimensions, hyperbolic 
spacetime with constant curvature —1 is given by the metric 

(4) ds2 = L^(da;2_^y2N 

where — 7r/2 < y < tt/2. Spherical spacetime is realized on the same manifold by 
simply inverting the sign in 

Consider a compact domain D G M . For the applications we have in mind it is 
sufficient to consider very "nice" domains, without global pathologies, or intricate 
smoothness issues, so we will assume for simplicity that the interior of D is an m 
dimensional, contractible, submanifold of M and that the boundary of D is C°°, 
except possibly at a lower dimensional C°° submanifold. The causal relation (past- 
future) is defined on pairs a,b G D hy a ^ b iS there is a future-pointing timelike 
curve in D from a to b. We will also assume that the relation ^ is a partial order 
relation, namely, there are no points a ^ b € D with future-pointing timelike curves 
(in D) from a to 6 and from b to a. 

Given a spacetime M and a pair of points a ^ b, a natural domain to consider 
is the interval [a, b] consisting, of all points z which satisfy a ^ z ^ b. The interval 
[a, b] for a pair of points in the Minkowski plane is shown in figure [U as the tilted 
square |a: -|- y| < S, whose vertices are a,b,e,f. The same domain doubles as 
the interval of [a, 6] in spherical spacetime if (5 < 7r/2 and of [e, /] in hyperbolic 
spacetime (restricting to \y\ < 7r/2). The interval [e,/] is a compact domain in 
hyperbolic spacetime if 5 < tt/2. 

A discrete spacetime associated with a compact domain D is a random partially 
ordered set on n points obtained by sampling the n points from D using the volume 
form and inducing the relation -< on the chosen points. 

We are interested in the asymptotic properties of discrete spacetimes associated 
with a domain D as n becomes large. In particular, we may think of a discrete 
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Figure 1. An interval [a,b] in the Minkowski plane. 

spacetime as a finite approximation to D. It is then natural to ask whether we can 
"reconstruct" D from larger and larger discrete spacetimes associated with it. 

At first this seems difficult since we have thrown away the original metric and 
retained only the causal relation and the volume form. However, results of Hawking, 
King and McCarthy, |28| . and of Malament, [38] show that, in the continuous 
setting, knowledge of causality and the volume form is essentially equivalent to 
knowledge of the metric. 

We can consider more closely the type of quantities that we may want to re- 
construct. Given a timelike curve C we may define its length L{C), better known 
as proper time, to be /^^ ds. In relativity theory, this parameterization independent 
quantity measures the time which passes on a clock which is attached to a point-like 
object whose trajectory is given by the curve C. Proper time is also the action for 
the trajectory of a point-like object between two points a,b £ M . This means that 
a free falling point-like object moving between events a and h will follow a critical 
path for this action, this is Einstein's equation of motion. The critical paths, which 
are maximizing in this case, are the timelike geodesies in analogy with the situation 
in Riemannian geometry (where they are minimizing). 

Assume that a,b € D, with a < b. We may define their "distance" in D to 
be LD{a,b) = sup,^ L{C), where the supremum is taken over all timelike curves 



4 



EITAN BACHMAT 



c G D, whose initial point is a and whose endpoint is b. It is known that in (a, b) 
is finite when D is compact and is achieved by a timehke curve, |41] . We shaU call 
such a curve a maximal curve. A maximal curve may have non-geodesic segments 
on the boundary of D, but in the interior of D it must satisfy the geodesic equation. 
It must also be differentiable at all points where the boundary is differentiable. 

We would like to have discrete analogues of the concepts of length and maximal 
curves. Given a partial order -< , a chain of size fc is a set of elements xi, . . . ,Xk such 
that Xi -< Xi^i for all 1 < i < fc. A chain in a discrete spacetime P associated with 
D is the obvious discrete analogue of a timelike curve. Under this correspondence 
the analogue of curve length is simply the size of the chain, and the analogue of a 
maximal curve between a and b will be a maximal sized chain xi, . . . ,Xk in P such 
that a ^ xi and Xk -< b. We let Lp{a, b) denote the length of such a maximal sized 
chain. 

We will say that an event occurs with high probability (w.h.p.) if it occurs with 
probability approaching 1 as n tends towards infinity. The following folk theorem 
provides the asymptotic high probability link between the discrete distance Lp{a, b) 
and the continuous distance L^^a, b). 

Theorem 1.1. Let D be a compact spacetime domain with the metric normal- 
ized so that Yol(D) = 1. Let a,b G D be such that a <b. Let P denote a discrete 
spacetime associated with D, then, for any e > 0, w.h.p. 



where Cd is a constant which depends only on the dimension. Moreover, given any 
Riemannian metric which is compatible with the differential structure of M , and 
any e > the maximal chain in P between a and b will be contained w.h.p. in an 
e -neighborhood of a maximal curve in D between a and b. 

The result has a somewhat tortured history. It is implicit in j40) . It was 
proved in Minkowski space in [17] , stated explicitly in [18j and proven in dimension 
d + 1 = 2 in [22j . The proof of [22) , which is in a completely different context (no 
mention of spacetime geometry), generalizes to arbitrary dimension. Meanwhile, in 
the physics literature on discrete gravity it seems to retain the status of a conjecture, 
[491 129) . A numerical study can be found in [29] . 

The constant Cd which appears in the theorem is known only for d = 1, where 
according to results of Kerov and Vershik and of Logan and Shcpp, |51L 133] we 
have ci = V2. For d > 1 there are only bounds for Cd, [17^ . which are obtained via 
first moment methods. 

The proof of theorem II. II is fairly elementary. It is first proved for an interval 
in Minkowski space. In this case it follows from a simple scaling argument and the 
transitivity of the Lorentz transformation group (the isometries) on timelike lines. 
The move from the local result to the global one is achieved using finite coverings 
of D by small patches where the metric is nearly constant, and by using the upper 
semi-continuity of the length functional, |41) . 

Given the law of large numbers Lp{a,b) ~ C£)i_D(a, &)^^^^'''^^'', it is natural 
to ask for a central limit theorem. Therefore, we define the random variable error 
term 



(5) 



Lp{a, b) - CdLD{a, b)n^/'^ 



(6) 



AD.n(a, b) = Lp{a, b) - coLoia, b)n^'^''+^^ 
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A general concentration result of Talagrand, }46|. shows that w.h.p. Var(AD_„) = 
(^(■yji/ld+i)-) Pqj. > 1^ essentially nothing is proven, beyond this fact. One expects 
that when the maximal curve is unique and entirely contained in the interior of D, 
the order of magnitude of l^D,n will only depend on the dimension of D. One may 
further expect that as n goes to infinity, l^D,n will have the limiting form 

(7) nP-'^'^+^hD 

where Pd is a dimension dependent constant and xd is some distribution. When 
d — 1, and D = [a, b] is an interval in the Minkowski plane, the famous Baik-Deift- 
Johansson theorem, [11] states that A_d_„ does have a limiting form as above with 
/3 = 1/3 and xd = F2 where F2 is the Tracy- Widom distribution of fluctuations of 
the largest eigenvalue in the GUE ensemble. Results of Baik and Rains \12l 113] 
show that Xd for the triangular domain a,e,b in figure [T] is F4, the Tracy- Widom 
distribution of GSE, and the distribution corresponding to the triangle, a, e, / is i^i, 
the Tracy- Widom distribution for GOE, see |47j for more on these distributions. 
These results are much deeper than the other results presented in this survey. 

Figure[5]shows the distribution of A^i^^/n^/^ in discrete spacetimes with 500, 000 
points associated with a hyperbolic, flat and spherical interval, which are repre- 
sented by the domain |a;| -I- \y\ < tt/S. Each distribution was composed of roughly 
350,000 samples. For the domain in question one can calculate CdLjj analytically 
to obtain CdLo ~ 1.973-\/n for the hyperbolic metric, CdLjj « 2.026y/n for the 
spherical and CdLu = 2 for the flat case. The figure also displays the average, 
variance, skew and kurtosis, for the distributions. The respective values for F2 (flat 
space asymptotics) are approximately —1.771,0.902,0.224,0.093, |47| . As can be 
seen, the numerics suggest that the following hypothesis: 

In hyperbolic and spherical spacetime the distribution xd has the form aF2+j, 
where 7 has the same sign as the curvature and a = 1 if the interval domain is 
calibrated so that CdLu = 2. A relation between the shift and curvature is expected 
since the curvature controls the length of curves which are close to the geodesic. 

Moving to d > 1 nothing is proven even for the flat interval case. The behavior 
of AD,n has been studied in the context of statistical physics because one of the 
simplest surface growth models, known as Polynuclear growth (PNG) can be de- 
scribed in terms of a discrete spacetime |43L I42j . In the PNG model one considers 
a crystal in contact with a super-saturated vapor. Once in a while a nucleus is 
formed on the existing surface and one assumes that the nucleus spreads evenly 
in all directions at a constant speed, forming a round new layer. When two such 
layers collide they coalesce. In terms of the PNG model, theorem 11.11 describes 
the macroscopic shape of the surface while the properties of Ajj^n describe finer 
aspects such as roughness of the surface. PNG is a model within a larger class of 
Kardar-Parisi-Zhang growth processes which are expected to exhibit some univer- 
sal features, including the values of the exponents Pd- Some heuristic theoretical 
considerations in physics, |32j , suggest that (32 and /Js have the form x/{2 ~ x) 
for some x of the form 2/(fc + 2), where k is an odd integer. Numerical measure- 
ments predict that (32 = 0.24, (i^ = 0.165, [42j . The closest values predicted by the 
heuristic argument would be /32 = 1/4 (fc = 3) and /^a = 1/6 (fc = 5), however, it 
is claimed in |42j that the numerical value (32 = 0.24 seems significant enough to 
exclude (32 = 1/4. A very interesting prediction of the physical heuristic arguments 
is that Var(A£i_„) < for all e > and n large enough, when d > 4 (essentially 
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Figure 2. A„ for intervals in flat, spherical and hyperbolic space. 

/3d = for d > 4.). The existence of such a critical dimension is subject to great 
debate in the physics literature, |20l 1321136] . and is currently not supported by nu- 
merical evidence. Numerics lead to an estimated value of P4 in the range 0.1 — 0.15. 
We note that many of the numerics are not computed on discrete spacetime, but 
rather, on closely related partially ordered sets in which the light cone is replaced 
by a "simplicial light cone" . The main example is the unit m dimensional cube with 
the partial order relation (xq, • ■ • , Xd) -< (xg, . . . , x'^) iff xt < x\ for all < z < d. It 
is expected that, after proper normalization, the fluctuations of the longest chain for 
this partial order are the same as those of the flat interval spacetime (universality) . 
Numerical calculations are easier and faster for the simplicial partial order. 

Another result which is strongly related to theorem 11.11 concerns the behavior 
of a random chain of length n, as n becomes large. We will state a semi discrete 
version. 

Theorem 1.2. Let a,b Cz D Let P = {xi, . . . Xn} be a discrete spacetime asso- 
ciated with D which is conditioned to be a chain between a and b , i.e., a ^ xi < 
X2 ^ ■ . ■ < Xn < b, then for any e > 0, w.h.p., P is contained in an e neighborhood 
of a maximal curve between a and b. 

Locally, where maximal curves and timelike geodesies coincide we can think 
of theorem 11.21 as providing a nice probabilistic interpretation for Einstein's law of 
motion, namely, a completely random trajectory of a particle between two endpoints 
is likely w.h.p. to closely follow a geodesic. So, perhaps God does play dice but as 
stated by a wise old man, "everything is foreseen, yet freedom of choice is given" , 
Rabbi Akiva, Avot 3:19. 
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Theorem 11.21 was proven in [221 , again in a very different context, for two 
dimensional domains, the proof extends to arbitrary dimension. A probabihstic 
interpretation along these lines to Einstein's equation of motion was sought in [24J . 
For similar and related results in the context of relativistic Markov processes and 
Brownian motion see [251 1261 127j . 

So far we have considered maximal chains in a discrete spacetime. We can 
consider the analogues problems for anti-chains. An anti-chain is a set of points no 
pair of which is comparable w.r.t. ^. In a continuous spacetime this is known as 
an achronal set. In the same way that chains are related to the concept of timelike 
curves, anti-chains are closely related to the concept of a spacelike hypersurface. 
A differentiable hypersurface is said to be spacelike if all its tangent vectors are 
spacelike. We have noted that geodesies locally maximize proper time. Proper 
time may be identified with the induced volume of a timelike curve. Similarly we 
may consider maximal hypersurfaces which are spacelike hypersurfaces which locally 
maximize the induced volume among spacelike hypersurfaces. The corresponding 
Euler-Lagrange equation which is satisfied by a maximal hypersurface states that 
the scalar curvature vanishes along the hypersurface. In analogy with theorems ll.il 
we can state the following result. 

Theorem 1.3. Given P tobe a discrete spacetime associated with a domain D 
of unit volume, let a{P) denote the size of a maximal anti-chain in P. Let a{D) 
denote the maximal volume of an achronal spacelike hypersurface in D, then there 
exists a constant such that for all e > 0, w.h.p., 

I Zd/idTT) 



n 



Moreover, if A d P is an anti-chain of size a[P), then, w.h.p., there is a spacelike 
hypersurface N of volume a{D) such that A is in an e neighborhood of N. 



Similarly, we have the analog of theorem 1 1.2) which states that a large random 
anti-chain in D will cluster around an achronal spacelike hypersurface of maximal 
volume. 

Geodesies and maximal hypersurfaces are examples of critical points of the 
volume functional. The Nambu-Goto action in string theory is also given by the 
volume functional, this time applied to surfaces which have a timelike tangent and 
a spacelike tangent at each point. Such surfaces describe the (classical) trajectory 
of a string. As before, one is interested in the critical points of this action. It would 
be interesting to provide some probabilistic interpretation to the critical points via 
causality preserving maps from a family of discrete spacetimes. 

We can use theorems 11.11 or 11.31 to provide a very simple procedure for recon- 
structing the dimension of D from a large associated discrete spacetime. Several 
solutions for this problem were proposed, [39L 1491 144] . however, all of them rely 
on specific properties of intervals in Minkowski space and do not work for general 
domains in general spacetimes. 

Let P be a given discrete spacetime associated with D. Randomly subdivide P 
into two disjoint sets Pi and P2 of size 2n/3 and n/3 respectively. We compute the 
sizes fci, k2 of longest chains in Pi, P2 and the sizes h, I2 of the maximal anti-chains. 
By theorem II. II the ratio fci/fc2 will approach 2^/^"^+^) w.h.p., therefore we can use 
1/ log2(A:i/fc2) as our estimate of the dimension of D. Similarly, by theorem 11.31 
the ratio h/h will approach 2'^/^'^+^) and so 1/(1 — 1/ log2(ii/?2)) can serve as a 
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dimension estimator. The basic idea of dividing P into two sets of uneven size is 
rather flexible and can be used in conjunction with many other properties of the 
partially ordered set in a similar way. The accuracy of the method depends on the 
properties of error terms like A^^n or a{P) — c'j^a{D)n'^/'^'^^^\ to the best of our 
knowledge the latter has not been explored even numerically. 

Having considered chains and anti-chains we may also consider a well known 
"bottom to top" peeling process which involves both. Let P be a finite partially 
ordered set. We denote by Ai the set of minimal elements of P, by A2 the set of 
minimal elements of P—Ai and more generally, by Ai the set of minimal elements of 
P — UAj , j < i. The number of such sets is well known to be equal to the size of the 
longest chain, hence we can analyze its asymptotic behavior using theorem ll.il The 
peeling process will figure prominently in the applications which will be described 
below. 

We can symmetrically define the "top to bottom" peeling process with A'^ 
defined recursively as the set of maximal elements. We call the set A[, consisting 
of the maximal elements in P, the maximal layer. The maximal layer plays an 
important role in computational geometry, especially in the two dimensional case. 
Consequently, the asymptotics of the maximal layer have been studied in depth in 
the case d = 1, [231 110] . We can consider D,nax, the set of maximal points of the 
domain D. Generically, D,nax will be a spacelike curve. We can consider its length 
L{Dmax) w.r.t. the metric —g for which it is timelike. We can restate the main 
result of |23j as saying that if vol{D) — 1, w.h.p. the size of the maximal layer 
will be approximately ^J^L{D„Lax)^/n, if L{Dmax) > 0. This result generalizes 
to any dimension if we replace L{Dmax)^y yo\{Dmax) and \/tt/2 by a constant, 
which is currently unknown. It is shown in [10] that the error term has a normal 
distribution with order of magnitude n^/^. The case where VoliLmax) — 0, as with 
intervals, has been studied in the context of the simplicial light cone partial order 
on the cube, [9l I15L I16j . The expected size of the maximal layer in this case is 
■jy^j^jy \n{n)"^~^ , and the fluctuations are again normal. We can expect similar 
results to hold for discrete spacetime, with different constants. The size of the 
maximal layer in an interval is of interest since it determines the number of edges 
in the graph representation of a discrete spacetime. 

Paper [23] also presents some examples of finite volume, non-compact domains, 
for which the above asymptotic estimate does not hold. From the point of view 
of spacetime geometry, the problem is that in the counter examples the length 
function w.r.t. —g is not upper semi-continuous at D^axj a problem which does 
not occur when D is compact. 

For additional reconstruction results the reader may consult |37j . 

2. A list of problems 

We present a few problems from different fields of applications. 

1) Scheduling of I/O (input/output) requests in disk drives: Modern disk 
drives such as the one in a PC or laptop can handle multiple requests to read and 
write data. They can rearrange the order of requests so as to handle them more 
efficiently. In practice this is a very important feature and in some cases can double 
the performance of the disk drive. The problem of finding the optimal order of the 
requests is known as the disk scheduling problem. We may also ask for the amount 
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of time that will be required to handle all the requests. This is known as service 
time. Since disk drives rotate at a constant speed we may ask equivalently for the 
number of disk rotations which are required to handle all the requests. We refer 
the reader to [2l 1481 [30j for more information on disk scheduling. 

As a toy model for this problem we may consider the simplistic case in which 
we assume that the read/write head of the disk accelerates and decelerates instan- 
taneously. Under these assumptions Andrews, Bender and Zhang (ABZ) gave an 
algorithm which provides an optimal ordering of requests, [2]. 

What is the expected service time of the ABZ algorithm for n requests? 

2) Connecting-the-dots via Lifschitz functions: The following problem was 
introduced in [3 in the context of estimating parameters for a multi-scale pattern 
recognition algorithm. 

Given a set of n points sampled via some density in a bounded domain in d+1 
dimensional Euclidean space. 

What is the maximal number of points which can lie on the graph of a Lifschitz 
function (on the first d coordinates) with Lifschitz constant 1? 

3) Airplane boarding: We consider passengers boarding an airplane. The input 
parameters for the process are 

• u - the average amount of aisle length occupied by a passenger; 

• w - the distance between successive rows; 

• b - the number of passengers per row; 

• D - the amount of time (delay) it takes passengers to clear the aisle, once 
they have arrived at their designated row; 

• An airline boarding policy. For example, the policy may first allow pas- 
sengers from the back third of the airplane to board, followed by the front 
third and finally the middle third. Another choice which is very popular 
is to board passengers from the back to the front, say, the back third, 
followed by the middle third and finally the front third. 

Given u, w, 6, D what is a good boarding policy for an airline to pursue? 

3. Problem solutions 

We show how the problems presented above are linked to discrete spacetime. 

1) Scheduling of I/O (input/output) requests in disk drives: Disk locations 
can be represented in polar coordinates (r, 9) , say with 1 < r < 2 and < < 27r 
and (r, 0) identified with (r, 27r). We are also given a distribution p{r,9) which 
reflects the popularity of different data in the disk. We can then assume that the 
data requests are sampled w.r.t. p. We assume for simplicity that the maximal 
radial speed of the disk is 1. We may then consider the Lorentzian metric 

(9) ds^ = 2p{r,0){dr^ -d9^) 

on the disk. Spacelike curves w.r.t. this metric correspond to possible motions of the 
disk head in a single rotation. Since the volume form for this metric is proportional 
to p(r, 0) we can think of n data requests as forming a discrete spacetime. We may 
apply the peeling process to this discrete spacetime. The ABZ algorithm has a very 
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natural interpretation in terms of the peeling process. The i'th layer At is precisely 
the set of requests that will be handled in the i'th rotation of the disk. The service 
time will then be the number of layers in the peeling process which according to the 
results of the previous section is asymptotic to L{C)y/n where L{C) is the length 
of the maximal curve in the model. Since the boundary of the model is spacelike, 
the maximal curve will be a geodesic in this case. In disk drives, files are laid out in 
consecutive tracks, each track corresponding to a certain fixed radius r. This type 
of data arrangement naturally leads to a popularity density p which depends only 
on r, namely p — p{r). In this case symmetry considerations show that the maximal 
curve is any curve which fixes the angle 9 and calculating the length we get that 
the service time concentrates around 'n/S/q^ \/p{r)dr^/n. Since we have an infinite 
number of maximal curves the magnitude of the error term Ao,n can be somewhat 
larger than expected. When the distribution p is the uniform distribution, it can be 
shown using results from |31|, I35j that An.n has order of magnitude n^/^ \og{n)^^^. 
We expect that the variance will actually be smaller in this case, in analogy with 
extreme value distributions. For more information, see [6l [7]. 

We need to say a few words about the (non) realism of this model. The as- 
sumption that the disk does not accelerate or decelerate is far from true. The 
model which we have used predicts that a disk can handle an order of magnitude 
of ^/n requests per rotation, a real disk will not handle more than 2-3 requests 
per rotation regardless of the number of requests, due to some fixed time delays. 
Nonetheless, the type of algorithms that are used in this case (ABZ) are indicative 
of successful algorithms in more complex and realistic models, therefore one does 
gain an understanding of the full problem by considering these simple cases. 

2) Connecting-the-dots via Lifschitz functions: The relation between 
this problem and spacetime geometry becomes clear once we note that by definition 
the graph of a Lifschitz function is a spacelike hypersurface in Minkowski space. 
In order to think of the points as points in a discrete spacetime we have to match 
the volume form of the metric to that of the point distribution so we introduce the 
conformally fiat metric 

(10) ds^ = p{x, t)2/('i+l) (di2 _ ^ ^^2y 

i 

By theorem 11.31 the largest number of points will be picked in an e neighborhood 
of a maximal hypersurface in the domain and the number of points will be pro- 
portional to the volume of the hypersurface. We also notice that this problem is a 
generalization of the previously discussed disk scheduling problem. 

3) Airplane boarding: We represent passengers by their position in the 
boarding queue x and their row number y, both normalized to produce a point 
(x, y) in the unit square. The airline boarding policy (and the way passengers react 
to it) leads to a distribution p{x,y) on the joint row/queue position. For example, 
assume that the policy is to board the back third of the airplane, followed by the 
front third and finally the middle third. The corresponding distribution p{x^ y) 
is given by p(x,y) = 3 on the sub-squares [0,1/3] x [2/3,1], [1/3,2/3] x [0,1/3], 
[2/3,1] X [1/3,2/3], and p{x,y) = elsewhere. The back-to-front poHcy with 3 
equal group sizes will lead to a distribution with p{x, y) = 3 on the subsquares 
[0,l/3x [0,1/3], [1/3,2/3] X [1/3,2/3] and [2/3,1] x [2/3,1]. 
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We define k = k measures the congestion of passengers in tlie airpiane. It 
is tfie reciprocal of the portion of aU passengers which can comfortably stand along 
the entire length of the airplane aisle. In real airplanes, k equals roughly 4. We 
also define a(x, y) = J p{x, z)dz. a measures the density of passengers which want 
to pass by the passenger {x,y) to get to their seats. Given p and k we define the 
Lorentzian metric 

(11) ds^ = AD'^p{x, y){dx(dx + ka{x, y) dy)) 

Since p{x, y) is proportional to the volume form, we can think of the passengers 
as a discrete spacetime sampled from the unit square equipped with this metric. 
It can be shown, [4] that the causal relation -< associated with the metric (fTTj) 
asymptotically coincides with the notion of passengers blocking other passengers 
from getting to their seats. Consequently, the boarding process roughly coincides 
with the bottom-to-top peeling process applied to this discrete spacetime. The 
first layer Ai , corresponds to passengers which can sit down without being blocked 
by other passengers. More generally, the set Ai asymptotically corresponds to the 
passengers which can sit unobstructed once the passengers in yjAj, j < i, have 
sat down. The boarding time is then given by L{C)y/n, where C is the longest 
(maximal) curve in the unit square with the metric (jlll In figures [3] and 21 we see a 
comparison of the maximal curve, with the longest chain of passengers blocking each 
other, for a huge airplane with 10^ passengers and congestion k = 0.5 and fc = 5 
respectively. We must admit that when the number of passengers is 200 or so the 
maximal curve does not follow the geodesic well. Furtheremore, it turns out that 
large portions of the maximal curve typically lie on the boundary of the support of 
p{x, y). This causes significantly greater variations in boarding time and makes the 
error term Ao,n much larger, [4]. Nonetheless, rather miraculously the ranking 
of boarding pohcies which is produced by the estimate £(C)i/n nearly coincides 
with the rankings given by very detailed computer simulations of the boarding 
process, |5l I50j . The methods above led to a detailed understanding of back-to- 
front policies as practiced by many airlines. Such policies are parameterized by the 
infinite dimensional simplex, consisting of sequences ro = < ri < . . . < Tm-i < 
1 = r„j, for any positive integer m. Given such a sequence, n passengers, with d 
passengers per row, the corresponding policy boards passengers in rows rm-in/d 
to Tmn/d = n/d, followed by passengers in rows rm-2'n/d to Vm-in/d and so on. 
Using these methods it was shown in [8] that for A; > 1 the expected boarding time 
of any back-to-front policy is at most 

faster than random boarding (no airline policy). For the realistic value of fc = 4 it 
shows that we cannot obtain a speed up of more than 20% in boarding time using 
back-to-front and in fact most back-to-front policies yield performance which is 
worse than random. The estimate [T^] which is based on constructing long (but not 
maximal) curves whose length satisfies a recursive inequality, is not tight. We are 
currently exploring more sophisticated recursions which will yield better bounds 
for any given k. We conjecture that for /c = 4 no back-to-front policy will yield a 
reduction in boarding time of more than 5% . For popular accounts of the airplane 
boarding problem see |14L 1341 152j . 
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Figure 3. The maximal curve and the longest chain for 10^ pas- 
sengers, when fc = 0.5 
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